Interplay of disorder and magnetic field in the superconducting vortex state 
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We calculate the density of states of an inhomogeneous superconductor in a magnetic field where 
the positions of vortices are distributed completely at random. We consider both the cases of s-wave 
and d-wave pairing. For both pairing symmetries either the presence of disorder or increasing the 
density of vortices enhances the low energy density of states. In the s-wave case the gap is filled and 
the density of states is a power law at low energies. In the d-wave case the density of states is finite 
at zero energy and it rises linearly at very low energies in the Dirac isotropic case {ao = t/Ao = 1, 
where t is the hopping integral and Ao is the amplitude of the order parameter) . For slightly higher 
energies the density of states crosses over to a quadratic behavior. As the Dirac anisotropy increases 
(as Ao decreases with respect to the hopping term) the linear region decreases in width. Neglecting 
this small region the density of states interpolates between quadratic and back to linear as ao 
increases. The low energy states are strongly peaked near the vortex cores. 

PACS numbers: 74.25.Qt, 74.72-h 



I. INTRODUCTION 

The interaction between the superconductor quasi- 
particles and the vortices induced by an external mag- 
netic field has been a subject of considerable recent 
interes^^. In the presence of vortices the quasiparticles 
feel the combined effect of the external magnetic field and 
of the spatially varying field of the chiral supercurrents. 
By performing a gauge transformation to effectively re- 
duce the system to the one in a zero average magnetic 
field it was shown^ that the natural low energy quasi- 
particle modes are Bloch waves rather than the Dirac 
Landau levels proposed in Refs. Q] and 0- The results 
or Ref£ also showed that the quasiparticles besides feel- 
ing a Doppler shift caused by the moving supercurrents^ 
(Volovik effect) also feel a quantum "Berry" like term due 
to a half-flux Aharonov-Bohm scattering of the quasipar- 
ticles by the vortices. 

The effect of disorder on the low-energy density of 
states of superconductors has also been a subject of much 
recent activity^, particularly in the case of <i-wave sym- 
metry. Several conflicting predictions have appeared in 
the literature which have mainly concentrated on the 
effect of the presence of impurities. Some progress to- 
ward understanding the disparity of theoretical results 
has been achieved realising that the details of the type of 
disorder affect significantly the density of states^. Par- 
ticularly in the case of <i-wave superconductors, in con- 
trast to conventional gapped s-wave superconductors, the 
presence of gapless nodes is expected to affect the trans- 
port properties. Using a field theoretic description and 
linearizing the spectrum around the four Dirac- like nodes 
it has been suggested that the system is critical. It 
was obtained that the density of states is of the type 
p(e) ~ |e| Q , where a is a non- universal exponent depen- 
dent on the disorder, and that the low energy modes are 
extended states (critical metal) 6 . Taking into account 



the effects of inter-nodal scattering (hard-scattering) it 
has been shown that an insulating state is obtained in- 
stead, where the density of states still vanishes at low 
energy but with an exponent a — I independent of 
disorde»i. The addition of time-reversal breaking cre- 
ates two new classes designated spin quantum Hall effect 
I and II, due to their similarities to the usual quantum 
Hall effect, corresponding to the hard and soft scattering 
cases, respectively^. The proposed formation of a pairing 
with a symmetry of the type d + id breaks time-inversion 
symmetry 8 but up to now remains a theoretical possibil- 
ity. On the other hand applying an external magnetic 
field naturally breaks time-reversal invariance and there- 
fore it is important to study the density of states in this 
case. 

In general, disorder is due to the presence of impuri- 
ties which may either scatter the quasiparticles and/or 
may serve as pinning centers for the field induced vor- 
tices. The density of states of a dirty but homogeneous 
s-wave superconductor in a high magnetic field, where 
the quasiparticles scatter off scalar impurities, was con- 
sidered using a Landau level basiA For small amounts 
of disorder it was found that p(e) ~ e 2 but when the 
disorder is higher than some critical value a finite den- 
sity of states is created at the Fermi surface. In the same 
regime of high magnetic fields, but with randomly pinned 
vortices and no impurities, the density of states at low 
energies increases significantly with respect to the lat- 
tice case suggesting a finite value at zero energy^. Refs. 
ITU and considered the effects of random and statisti- 
cally independent scalar and vector potentials on e?-wave 
quasiparticles and it was predicted^ that at low energies 
p{e) ~ po + ae 2 , where po ~ B 1 / 2 . The effect of randomly 
pinned discrete vortices on the spectrum of a d-wave su- 
perconductor was considered recently and a preliminary 
report was presented in Ref. for the isotropic case. 
In this work we study the density of states as a function 
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of vortex density and consider the effects of the Dirac 
anisotropy. We also study the local density of states 
(LDOS) and the inverse participation ratio (IPR) and 
study the effect of disorder on the localization of the low 
energy states. 

The nature of the low energy states in the presence of 
a single vortex with s-wave symmetry was solved long 
agoii There are localized bound states in the vortex 
core. The d-wave symmetry case led to some early con- 
troversy but it was eventually clearly demonstrated that 
the low energy states, even though strongly peaked near 
the vortex core, extend along the four nodal directions 
and are indeed delocalized as shown by the behavior of 
IPRi£. In the vortex lattice the states are naturally also 
extended^*i&. In the clean limit it is therefore expected 
that the external field will increase the low-energy den- 
sity of states. The expectation that these are delocalized 
is evidenced by the increase of the thermal conductivity 
with fielcUi in contrast to the reducing effect of conven- 
tional superconductors^. 

The addition of impurities in zero magnetic field has 
been studied using the Bogoliubov-de Gennes (BdG) 
equations. It was found that the d-wave superconduc- 
tivity is mainly destroyed locally near a strong scat- 
terer. The superfluid density is strongly suppressed near 
the impurities but only mildly affected elsewhere^. No 
evidence for localization of the low energy states was 
found and accordingly the superfluid density is indeed 
suppressed but less than expected2£L2i and, accordingly, 
the decrease of the critical temperature with disorder 
is much slower than previously expected in accordance 
with experiments^. Similar results of an inhomoge- 
neous order parameter were also obtained for s-wave 
superconductors^. 

The question we address in this paper is the influence 
of the positional disorder of the vortices on the quasipar- 
ticle states of either a s- or a d-wave superconductor in 
an external magnetic field and we will not consider the 
scattering off impurities. We obtain that the low energy 
states are strongly peaked near the vortex locations as 
evidenced by the LDOS. At higher energies the states 
have a very uniform distribution throughout the system. 
For these states the IPR scales as 1/L 2 indicative of ex- 
tended states. In the case of the low energy states the 
same quantity does not follow this scaling but does not 
saturate either. This is probably a consequence of finite 
size effects indicative that the system sizes considered are 
smaller than the localization length. The results indicate 
at best a large value for the localization length but are 
more consistent with extended states rather than with 
localized states. 

In section II we describe the model and present the 
BdG equations to be solved and consider the effect of 
the vortices on the supercurrent profiles. In section III 
we study the effects of positional disorder on s-wave su- 
perconductors. In section IV we consider d-wave super- 
conductors for the isotropic and anisotropic cases. Also 
we study the spatial structure of the low-lying quasipar- 



ticle states. We end with the conclusions. 



II. DESCRIPTION OF THE MODEL 

The strong correlations of the high T c materials are 
typically modelled using a Hubbard-like Hamiltonian for 
the electrons. The superconductivity is considered using 
a BCS-likc formulation which provides good agreement 
with experimental results. Even though the normal phase 
of these materials is particularly challenging it is by now 
accepted that the phenomenology in the superconducting 
phase is reasonably well described by BCS theory. 

A. Bogoliubov-de Gennes Hamiltonian 

We will consider the lattice formulation of a disordered 
d-wave superconductor in a magnetic field. We start from 
the Bogoliubov-de Gennes equations Hip = ei\> where 
?/>t(Y) = (it*(r), v*(r)) and where the matrix Hamiltonian 
is given by 

*-(a.-*0 « 

with£J£ 

h = -tJ2e-^^ +SA{r) - dl ss~e F (2) 

8 

and 

A = A £V*M^ e **M < (3) 
s 

The sums are over nearest neighbors (8 = ±x, ±y on the 
square lattice); A(r) is the vector potential associated 
with the uniform external magnetic field B = V x A, 
the operator s$ is defined through its action on space 
dependent functions, ssu(r) = u(r + 8), and the operator 
77,5 describes the symmetry of the order parameter. 

The application of the uniform external magnetic field 
B generates in type II superconductors compensating 
vortices each carying one half of the magnetic quantum 

flux, — . We take the London limit, which is valid for low 

magnetic field and over most of the H — T phase diagram 
in extreme type-II superconductors like the cuprates, as- 
suming that the size of the vortex cores is negligible and 
placing each vortex core at the center of a plaquette (unit 
cell). The vortices are distributed randomly over the 
L x L plaquettes. In this limit we can take the order 
parameter amplitude Ao constant everywhere in space 
and we can factorize the phase of the order parameter as 
shown in Eq. |3J 

At this stage, it is convenient to perform a singular 
gauge transformation to eliminate the phase of the off- 
diagonal term J2J) in the matrix Hamiltonian. We con- 
sider the unitary FT gauge transformation H — > U~ l TW , 
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FIG. 1: Color vector field plot showing the profile of the 
conventional superfluid wave vector k s = ~ (kf + kf ) for a 
regular vortex lattice and for B = 1/144. 



FIG. 2: Color vector field plot showing the profile of the 
conventional superfluid wave vector k 3 = | (kf + kf ) for a 
random configuration of vortices and for B = 1/144. 



where' 

with 4>a(?) + 4>b(y) — 4>( r )- The phase field <p(r) is 
then decomposed at each site of the two-dimensional 
lattice in two components 4>a(y) and 0b (r) which are 
assigned respectively to a set of vortices A, positioned 
at {rf}i=i.N A > an d a set of vortices B, positioned at 
{rf }i=i,jv B - The phase fields §h=a.b are naturally de- 
fined through the equation 

V x V</y(r) = 27rz^()(r-if) (5) 

where the sum runs only over the /x-type vortices. After 
carying out the gauge transformation Q the Hamilto- 
nian (JTJ reads 

H ' = \ AoJ^e-^fjU^ tJ2e- iVi(r) §s+e F ' 

V 5 5 / 

(6) 

The phase factors are given byi 6 V^(r) = f* + ■ dl 
and <5</>(r) = 4>a (r) — 0b (r), where ftk^ = mv^ 1 = 
tiVcf)^ — -A is the superfluid momentum vector for the 
/i-supercurrent. Physically, the vortices A are only vis- 
ible to the particles and the vortices B are only visible 
to the holes. Each resulting /i-subsystem is then in an 
effective magnetic field 

Kn = - ™ V x = B - o z ]T 5 2 (r - rf ) (7) 

i 

where each vortex carries now an effective quantum mag- 
netic flux (f>o- For the case of a regular vortex lattic e 3 ' 16 , 



these effective magnetic fields vanish simultaneously on 
average if the magnetic unit cell contains two vortices, 
one of each type. More generally, in the absence of spa- 
tial symmetries, as it is the case for disordered systems, 
these effective magnetic fields 'B^ A:B vanish if the num- 
bers of vortices of the two types A and B are equal, i.e. 
Na — Ng , and their sum equals to the number of elemen- 
tary quantum fluxes of the external magnetic field pene- 
trating the system. As the number of vortices in the 
two-dimensional system of size L x L is proportional to 
the quantized magnetic flux piercing through the system, 
we choose to parametrize the magnetic field intensity by 
the ratio of the number of vortices, — Na + Nb, by 
the number of lattice sites, B = jp^r- 

We consider here the disorder induced by the ran- 
dom pinning of the vortices over the two-dimensional 
L x L lattice. As the effective magnetic fields {7J) expe- 
rienced by the particles and the holes vanish on average 
within the gauge transformation Q we are allowed to 
use periodic boundary conditions on the square lattice 
(ty(x + nL, y) = ^(x, y + mL) = ^(x, y) with n, m G Z). 
The L x L original lattice becomes then a magnetic su- 
percell where the vortices are placed at random with a 
mean intervortex spacing I — 1/yB. The disorder in the 
system is then established over a length L. 

In order to compute the eigenvalues and eigenvectors 
of the Hamiltonian © we seek for eigensolutions in the 
Bloch form *„ k (r) = e~ 4k ' r ([/* k , V* u ) where k is a point 
of the Brillouin zone. In the following we will label these 
eigensolutions with the index n = (n, k). We diagonalize 
then the Hamiltonian e~ ikr 7iV kr for a large number 
of points k in the Brillouin zone and for many different 
realizations (around 100) of the random vortex pinning. 
The results of these computations are shown in section 
IIIII for the s-wave disordered superconductor case and in 





FIG. 3: Quasiparticle density of states for the s-wave sym- 
metry and for a regular lattice of vortices (no disorder). The 
magnetic field is B = 1/18, the magnetic unit cell contains 
2 vortices piercing an area of 6 x 6. The parameters are 
|U = -2.2t and A = t. 



FIG. 4: Quasiparticle density of states for the s-wave sym- 
metry and for a disordered lattice of vortices. The magnetic 
field is B — 1/18 and the linear system size is L — 18. The 
parameters are fi = — 2.2t and Aq = t. 



scction llVl for the d-wavc disordered superconductor case. 



B. Profiles of supercurrents 

The /i-superfluid wave vector k£(r) characterizes the 
supercurrents induced by the /i-vortices. This vector can 
be calculated for an arbitrary configuration of vorticesi& 
like 



k?(r) = 2nj 



a Ik ikxz 



^2 e ik (r " r 



(8) 



Here A is the magnetic penetration length and the sum 
extends over the infinite number of /i-type vortex posi- 
tions. We consider the case A — * oo which is an excel- 
lent approximation in extreme type-11 systems like high 
temperature superconductors. In this limit the repulsive 
interaction between vortices is not screened and there- 
fore the vortex distribution is not strictly arbitrary - 
long range interactions will try to force the vortex system 
to take only incompressible configurations. For the pur- 
poses of this paper, we assume that the pinning centers 
are strong enough to overcome the vortex repulsion over 
the lengthscales relevant to experiments. Thus, we are 
considering the limit of strong pinning. We have checked 
that the case with a random distribution of vortices is 
qualitatively the same as for the case where the vortex 
positions are allowed to vary with a radius of a few unit 
cells around a regular lattice position. 

The /i-supcrfluid wave vector distribution is completely 
determined by the configuration of the /i-vortices, as can 
be seen from Eq. [S] and is independent of the pairing 
symmetry. Since we are taking the London limit, where 



the vortex core size is negligible, we neglect any possible 
different symmetry contributions from inside the vortex 
core" and the chiral supercurrents simply reflect the cir- 
culation of the superfluid density around the vortex. 

In order to compute efficiently the /x-superfluid wave 
vector we use the translational symmetry introduced by 
the periodic repetition of the L x L supercells and then 
we are able to use the Fourier series representation of 
expression ©, e.g. 



2iir x -> 



N u 



Q X Z V c iQ-(r-rt ) 



Q 2 



E 

i=l 



(9) 



where Q = |? ("4, ^2) with n\^n<z € Z. We remark that 
the sum over vortex positions runs now only over the /Li- 
type vortices within a supcrccll and not, as in Eq. [5J over 
the infinite number of /x-type vortices present in the two- 
dimensional system. The price to pay for this procedure 
is the preceding infinite sum over the reciprocal vectors 
Q. We truncate this latter sum when the convergence 
of each components of k^(r) in each lattice point r is 
attained. 

Figures and show examples of distributions of the 
conventional superfluid wave vector which is half the sum 
of the two types of superfluid wave vector, 



k s (r) 



k»+kf(r) 



(10) 



In Fig. ^ we show the profile of the supercurrent veloc- 
ities in the lattice case. The unit cell has two vortices, 
one of each type A and B. In Figure 2 we show the su- 
percurrents for an arbitrary configuration of the vortices. 
For example, the 24 x 24 lattice shown in Fig. |2Jcan illus- 
trate a disordered supercell corresponding to a magnetic 
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FIG. 5: Quasiparticle density of states (s-wave) for different 
magnetic fields B = 2/200 (A), B = 4/200 (■), B = 7/200 
(□), B = 11/200 (•) andB = 20/200 (o) in units of hc/(2e8 2 ). 
The linear system size is L = 20 and the parameters are /i = 
— 2.2t, Ao = t. For clarity the different curves are vertically 
shifted. 

field B = 1/ 144 with two A- type vortices and two -B-type 
vortices. As it should be, the labelling of the A- and the 
-B-type vortices is completely arbitrary. The intensity of 
the supervelocities for the regular and the disordered vor- 
tex lattice peak around each vortex and decrease rapidly 
outside the vicinity of the vortex core (~ 55). 



III. DISORDERED s-WAVE 
SUPERCONDUCTORS 



or not, the tails of the quasiparticles spectrum (e 3> Ao) 
are expected to be described by quantized Landau levels. 

In Fig. 3 we show the density of states for a field 
B = 1/18 in the case of no disorder. The gap is clearly 
seen at low energies characteristic of s-wave symmetry. 
The states inside the gap (for |e| < t in Fig. |3J are the 
typical Caroli - de Gennes - Matricon (CdGM) bound 
states^. At higher energies the Landau levels are clearly 
visible. 

In Fig. 4 we show the effects of disorder for the same 
B field in a lattice of 18 x 18 sites (18 vortices). The 
gap is filled by the disorder. The reader should observe 
here that our disorder is in certain sense "infinite" since 
we place vortex positions completely at random. Thus, 
it appears that such full randomness in vortex positions, 
generated by strong pinning, suffices to close the gap in 
the single particle density of states. In constrast, the 
Landau level quantization structure clearly persists at 
high excitation energies. In general the increase of the 
magnetic field modifies the curvature of the quasiparticle 
density of states for the low energies. As it is shown in 
Fig. [S] the density of states seems to be decribed by the 
power-law formula 

p(e) ~ e a . (12) 

In Fig. El we fit the magnetic field dependence of the 
exponent a. This exponent obeys the following law 

a~c£-d (13) 

where £ = \j\fB is the mean intervortex spacing and d 
is found close to 1. In the cases of a strong magnetic 
field the density of vortices is high and strictly we are 
in a regime where the size of the vortex cores can not 
be neglected. In this regime the Landau level structure 
at high energies is clear and it extends to low energies 
superimposed by the effects of disorder. The lower limit 
for this high magnetic field regime is the value B ~ 0.16 
for which a — in Fig. [B] 



For the convcntionnal s-wave case the operator charac- 
terizing the symmetry of the order parameter is constant 
fjs = \ and the off-diagonal terms of the Hamiltonian © 
are then considerably simplified 

f - t J2^ (r) ss-e F Ao \ 

V 8 J 

(11) 

The s-wave Hamiltonian l|ll|) describes coupled parti- 
cles and holes evolving each in an effective magnetic field 
composed by the external magnetic field and the counter- 
acting field of cj>o flux carying vortices. The effective mag- 
netic fields are characterized by the Peierls phase factors 
y£~ A ' B . Independently of the system being disordered 



IV. DISORDERED d-WAVE 
SUPERCONDUCTORS 

For the unconventional d-wave case the operator fjg 
takes the form fjs = (—l) Sy is, we recall that §g acts on 
spatial dependent functions as §gu(r) = u(r+S) and that 
S = ±x, ±y characterizes unit displacements (hops) on 
the lattice. With these definitions the cf-wave Hamilto- 
nian can be derived from the Hamiltonian JSJ and reads 

/ -tJ2e^ {r) sg-e F A J2e lAs{r)+z ^sg\ 

V s s J 

(14) 




FIG. 6: Magnetic field dependence of the exponent a ex- 
tracted from the power-law fit p(e) ~ e a of the data presented 
in Fig. The symbols (A), (■), (□), (•) and (o) denote the 
same values of B as in Fig|S| the symbols (+) denote others 
values of B not presented in Fig|S] 

where the phase factor .4,5 (r) has the form 

Mr) = \[ +S - V0 B ) • dl 

= \[ +S (kf-kf)-dl. (15) 
In the Hamiltonian (|14|l and in Eq. ^|the vector 

a, = ±(k?-kf) (16) 

acts as an internal gauge field independent of the external 
magnetic field—. The associated internal magnetic field 
b = V x a s consists of opposite A — B spikes fluxes 
carying each one half of the magnetic quantum flux </>o, 
centered in the vortex cores and vanishing on average 
since the numbers of A- and B-type vortices are the same. 

A. Differences from the regular vortex lattice case 

The situation where the vortices are regularly dis- 
tributed in a lattice was treated before^ 6 ". Since the aver- 
age effective magnetic field vanishes it is possible to solve 
the BdG equations using a standard Bloch basis- the su- 
percurrent velocities are periodic in space and there is 
no need to consider the magnetic Brillouin zone. Taking 
the continuum limit and linearizing the spectrum around 
each node effectively decouples the nodes. It was shown 
that the low-energy quasiparticles are then naturally de- 
scribed as Bloch waves^ and not Dirac-Landau levels as 
previously proposed-' 2 . However, it was found that in the 
linearized problem different assignments of the A and B 




FIG. 7: Quasiparticle density of states for the d-wave case 
without (dashed line) and with disorder (solid line). The 
intensity of the magnetic field is B = 1/121, A = 0.5t, 
fi — —2.2t and the linear system size is L = 22. 

vortices lead to somewhat different spectra, which was 
unexpected^. It was found that defining theory on the 
lattice regularized this problem and indeed the system 
has a manifest internal gauge symmetry such that the 
spectrum is independent of the A-B vortex assignments, 
as it should be. Moreover, the lattice formulation ex- 
plicitly involves internodal contributions which are im- 
portant for the properties of the density of states in the 
disordered case. In the vortex lattice case, however, it 
was found that only in special commensurate cases (for 
the square lattice) the inclusion of the internodal contri- 
butions is relevant since only in such cases a gap develops 
due to the interference terms between the various nodes, 
estimated to be of the order of \[B. In a general incom- 
mensurate case the interference is not relevant leading 
to qualitatively similar spectra. In the d-wave case the 
spectrum is gapless with a linear density of states at low 
energy 16 . One would therefore expect that in a general 
disordered vortex case internodal scattering might not 
be relevant (particularly for high Dirac cone anisotropy 
an = vf/va = t/A > 1). 

In Fig. 13 we compare the densities of states for the 
lattice case and a case with disorder. At weak fields the 
density of states is small at low energies having a dip close 
to zero energy. We have checked for finite size effects on 
the spectrum. For system sizes larger than 16 x 16 the 
density of states at not very low energies converges and 
the finite size dependence is negligible. 



B. Isotropic case 

In this subsection we will focus our attention on the 
isotropic case qd = 1 in order to extract the general 
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FIG. 8: Quasiparticle density of states for the d-wave case 
and for different magnetic fields B = 1/200 (A), B = 3/200 
(A), B = 5/200 (■), B = 7/200 (□), B = 9/200 (•), 
B = 11/200 (o), B = 20/200 (♦) and B = 25/200 (0) in 
units of hc/(2eS 2 ). The linear system size is L = 20 and the 
parameters are p = — 2.2t, Ao = t. The inset shows the fits of 
the density of states (solid lines) inside the presented energy 
interval. For clarity not all the fits are shown. 



FIG. 9: Low-energy quasiparticle density of states p(e) for the 
d-wave case and for different magnetic field B = 1/200 (a), 
B = 3/200 (A), B = 5/200 (■), B = 7/200 (□), B = 9/200 
(•), B = 11/200 (o), B = 20/200 (♦) and B = 25/200 (0) in 
unit of hc/(2eS 2 ). For each data sets the linear system size 
is L — 20, the same as in Fig. |H| In the left panel the solid 
lines are linear fits of the dip region below e = 0.02f of the 
type p(e) — podip + /3\e\. In the right panel we present the 
near scaling at low energies (see text). 



dependencies of the quasiparticle density of states in the 
magnetic field B. 

In Fig. [5] we show the density of states for a system 
with size 20 x 20 and for various magnetic fields. The den- 
sity of states at small energies is finite up to quite small 
energies where there is a dip to a value that decreases as 
the magnetic field decreases. Only for quite small mag- 
netic fields the density of states approaches zero at the 
origin. Neglecting the narrow region close to the origin 
we have fitted the density of states using the power law 

p(e)=p + (3\e\ a . (17) 

In the inset of Fig. [S]we show the fits for the various val- 
ues of the magnetic field. Reasonable fits are obtained 
taking a ~ 2 and we obtain that po ~ B x l 2 . The various 
system sizes fit in the same universal curve indicating 
that the finite size effects are negligible. Note that in 
the lattice case the density of states at low energies is 
linea»iii£i2£iSi (this result differs from the behavior ob- 
tained by others for a <i-wave superconductor with no 
disordes^Si, where p(e ~ 0) ~ B 1 ^ 2 ). The finite den- 
sity of states at zero energy is therefore a consequence of 
finite disorder. 

In Fig. |2| we focus on the narrow region close to e = 
for the same set of parameters considered in Fig. |SJ Ex- 
cept for the lowest field case B = 1/200 the density of 
states seems to be finite at zero energy. Here the field 
density B — 1/200 corresponds to the particular case 
where only two vortices pierced the 20 x 20 disordered 
supercell. As shown in Ref. |25| in this case the spectrum 



is usually gapped and therefore the density of states van- 
ishes at zero energy. Performing a fit like in Eq. 1171 we 
obtain an exponent which is now close to 1. In this regime 
Podip also scales with y/~B and the slope scales linearly 
with B. In this low energy regime the finite size effects 
are still noticeable but the dependence on the magnetic 
field is common to the various system sizes. At these low 
energies the density of states for the various system sizes 
appears to be of the following approximate form 

p(e)~— tW— t), (18) 

where loh ~ V AB and J 7 is a universal function. In 
the left panel of Fig. [5| we show p(e) for various fields 
while in the right panel we illustrate the near scaling at 
low energies consistent with J-(x) ~ c\ + c-ix at small 
x. Note that f3(B — > 0) appears to be small but fi- 
nite, consistent with a crossover to a "Dirac node" scaling 
p(e) - (l/w ff )(l/P)J r (e/w H ) at very low fields. 



C. Anisotropic case 

We investigate now the dependence of the quasiparti- 
cle density of state on the value of the Dirac anisotropy 
ratio an- Such a study is interesting in order to compare 
our results with experiments; indeed for high-T c super- 
conductors such as YBa2Cu307 and I^S^CaC^Os the 




FIG. 10: Quasiparticle density of states for the d-wave case 
for different values of the anisotropy parameter: an — 1 (°), 
a D = 2 (•), a D = 3 (□), a D = 4 (■) and a D = 5 (V). The 
magnetic field intensity is B — 1/100. The linear system size 
is L — 20 and the chemical potential is p = — 2.2t. 



FIG. 12: Inverse participation ratio £ as a function of the 
quasiparticle energy e for B = 1/25 and for different linear 
system sizes L — 10 (o), L = 20 (♦) and L = 30 (A). Inset: 
scaling of the inverse participation ratio with 1/L for states 
in the close vicinity of the Fermi surface |e| < 0.05t (£o, B) 
and for states with energy |e| ~ t (£i.o, •)• 




=A ( /t 



FIG. 11: Fits of the data presented in Fig. ^] Upper left 
pannel: quasiparticle density of states presented in Fig. 1101 
(the same symbols are used) as a function of the rescaled 
energy e/Ao. Solid lines are fits of the power law type p(e) = 
po + F(e/Ao) a . On the others pannel solid lines are linear 
fits. Upper right pannel: zero energy quasiparticle density of 
states po as a function of \faE. Lower left pannel: Power 
law exponent a as a function of 1/ -J 'aD- Lower right pannel: 
Parameter T of the power law as a function of l/an- 



Dirac anisotropy ratio isiL2£ respectively an — 14 and 
old — 19. 

In the lattice case a high anisotropy increases the den- 
sity of states at low energies and leads to lines of quasi- 
nodesiSi2^. At high anisotropy (Ao « t) the nodes are 
very narrow and a one-dimensional like character is evi- 
denced. 

As shown in Fig. I1UI the low-energy quasiparticle den- 
sity of states at constant field is filled when the anisotropy 
parameter aD is increased. There is a narrow linear re- 
gion close to the origin that decreases as ac increases. 
Fig. ^2 shows the power law fits of the data presented in 
Fig. E]]neglectmg a very narrow region of width ~ 0.025£ 
around e = 0. It turns out that the low-energy quasipar- 
ticle density of states has the form 



p(e) ~ Po + r 



c 

Ao" 



with zero-energy density of states 



the exponent 





and the parameter 



Ao 

t 



(19) 



(20) 



(21) 



(22) 
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FIG. 13: Upper left panel: Color density plot of the quasiparticle local density of states p(r, e) for e ~ 0. Lower left panel: 
Color density plot of the quasiparticle local density of states p(r, e) for e ~ t. Upper right panel: Color vector field plot showing 
the profile of the internal gauge field a s = | (k^ — kf ). Lower right panel: Color vector field plot showing the profile of the 
conventional superfluid wave vector k s = | fk^ + kf). All the panels correspond to the same configuration of the vortex 
pinning and to the same magnetic field B — 1/144. 



The dependence of the exponent is interesting since in 
the isotropic case we retrieve a = 2 and with increasing 
anisotropy the exponent decreases to 1 (see lower left 
panel of Fig. Illjl characteristic of the Dirac nodes. At e ~ 
and for a high Dirac anisotropy ratio the quasiparticle 
density of states flattens since the coefficient T in <|19|) 
decreases with A . 



D. Spatial structure of the low-lying quasiparticle 
states 

As discussed above, the nature of the low lying states is 
an important point which is relevant due to the presence 
of disorder. A standard way to analyse the nature of 
the states is to study the IPR. The 1PR is defined in the 



usual way 

(E», r (l«n (r)| 4 + K (r)| 4 W-£ n )\ 

£(e) = - — " 1- (23) 

(En, r (K(r)| 2 + K (v)\ 2 )s(e-E n )) 

The brackets denote the averaging over disorder config- 
urations. The IPR £ is a direct measure of the spatial 
extend of the quasiparticle wavefunctions. It scales as 
l/L 2 for extended states and is constant for localized 
states with localization length £ c < L. 

Fig. E| presents the IPR £ as a function of the quasi- 
particle energy e and for different system sizes L. All the 
states are well extended, since there is no size dependence 
for the quantity L 2 £(e), except those in the close vicin- 
ity of the Fermi surface (|e| < 0.5i). In the inset of Fig. 
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FIG. 14: Quasiparticle local density of states p(r, e) for the 
d-wave case at different lattice sites r = (11, 11), (12, 12), 
(13,13), (14,14), (15,15), (16,16) around a vortex position 
r v = (13.5, 13.5) belonging to a regular vortex lattice. The 
magnetic field is B — 1/144. The other parameters are Ao = t 
and ii — —2.2t. 

El we present the scaling of the IPR £ with 1/L 2 . We 
consider an average over states at the Fermi surface with 
|e| < 0.05i (£o in the inset) and an average over states 
with e ~ t (£i.o in the inset). The states close to e = t are 
clearly extended since the scaling with 1/L 2 is quite ac- 
curate. The nature of the states close to e = is however 
less clear. From Fig. ^]we see that the quasiparticle 
states close to e = deviate from the strict 1 /L 2 scaling 
law (see the weak intercept of £o in the inset of Fig. I12fl . 
however our results show an obvious size dependence for 
£o and do not show a saturation of £o with the linear 
system size L. This latter fact indicates, in the hypoth- 
esis of an eventual localization of these low-lying states, 
that we are still far from the thermodynamic limit with 
the linear system sizes we can currently attain within our 
model. 

To gain further insight into the nature of the low lying 
states we show in Fig. ED the LDOS for the states close 
to e = and those close to e = t. At low energies the 
LDOS defined by 

p(r, e) = £ (Mr)| 2 + K(r)| 2 ) S(e - E n ) (24) 

n 

is strongly peaked near the vortex cores (upper left panel 
of Fig. U3J). At higher energies (e ~ 1) the LDOS 
is much more homogeneously spread over the system 
indicative of extended states; the values of the LDOS 
over the lattice are approximatively constant and slightly 
fluctuate around the expected value for extended states 
1/L 2 ~ 0.0017 (lower left panel of Fig. U3J). 

The right side panels of Fig. ^] show the profiles of 
the internal gauge field a s (upper panel) and of the con- 
ventional superfluid wave vector k s (lower panel) corre- 
sponding to the same configuration of vortex pinning as 



FIG. 15: Quasiparticle local density of states p(r, e) for 
the d-wave case at different lattice sites r = (3, 14), (4, 15), 
(5, 16), (6, 17), (7, 18), (8, 19) around a vortex position r v = 
(5.5, 16.5) belonging to a disordered vortex lattice. The mag- 
netic field is B = 1/144 and we use the 24 x 24 disordered 
supercell used in Fig. 1131 The other parameters are Ao = t 
and n = —2.2t. 

used in the computation of the LDOS (left side panels). 
For the profile of the superfluid wave vector k s the inten- 
sity of the supercurrents is peaked around each vortices 
and rapidly decreases as the inverse of the square of the 
distance to the vortices. The profile of the internal gauge 
field a s is more interesting. It depicts the two types of 
vortices A and B whirling in opposite directions (see up- 
per right panel of Fig. ILljl . Due to this opposite cir- 
culation around A- and B-type vortices, the interaction 
between a A- and a -B-type vortex produces significant 
field currents between the vortices. Fig. ^| shows a spa- 
tial correlation between these intervortex currents and 
the inhomogeneity of the local density of states between 
the vortices for the low-lying states e ~ (see upper left 
panel of Fig. I13fl . 

If we compare now the results for the LDOS with those 
for the IPR £ for the low-lying states, we can argue the 
following: as the participation ratio ^q 1 physically counts 
the number of lattice sites occupied by the quasiparti- 
cle wave function ^(r), and as the low energy states 
are mainly located around the vortices, the participa- 
tion ratio ^q 1 should scale for a fixed magnetic field as 
the number of vortices in a supercell ^q 1 ~ ~ BL 2 . 
This argument explain the fact that £o in the inset of Fig. 
1 121 seems to follow the 1/L 2 without saturation, the weak 
intercept being then negligible. The low-lying quasiparti- 
cle states appear then to be delocalized although strongly 
peaked around the vortex cores. 

We present now a spatial scanning of the LDOS in the 
vicinity of a vortex core for the case of a regular vortex 
lattice (Fig. I14f> and for the case of a disordered vor- 
tex lattice (Fig. I15fl . For both cases the LDOS away 
from the vortex core are qualitatively the same and are 
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comparable to that of a d-wave superconductor in a zero- 
magnetic field. Also for both cases the low-lying states 
are predominant in the close vicinity of the vortex core 
but their respective spectra are different. For the reg- 
ular vortex lattice case (Fig. I14f> we remark a double 
peak structure around the vortex core. This result is 
in qualitative agreement with the double peak in con- 
ductance observed in YBa2Cu307_,5 by scanning tunnel- 
ing microscopy (STM)2i. For the disordered vortex case 
(Fig. IT5jl zero energy peaks (ZEP) appear in the close 
vicinity of the vortex core (closest neighbor sites) and 
rapidily disappear (typically over 3(5) when moving away 
from the vortex core. We note also that close to the 
vortex cores the coherence peaks dissapear in both the 
regular and the disordered cases. 

V. CONCLUSION 

In summary, we have calculated the density of states 
of a disordered superconductor in a pinned fully random 
vortex array. Both the disorder and the magnetic field 
fill the density of states at low energies. In the s-wave 
case the density of states behaves as a power law with 
an exponent that scales with 1/yB. In general we find a 
finite density of states at zero energy for the d-wave case 
except in the limit of very small magnetic fields. The 
density of states deviates from the zero energy value by 
a power law. The zero energy density of states scales with 
the inverse of the magnetic length (VB). In the <i-wave 



case the Dirac anisotropy further increases the weight of 
the density of states at low energies. Also it affects the 
exponent of the power law. In the isotropic case the ex- 
ponent is 1 at very low energies and around 2 neglecting 
this narrow region. In the linear regime the density of 
states is of the form of a scaling function of the energy 
and the magnetic field. As the anisotropy increases this 
narrow regime shrinks considerably and, neglecting this 
region, the exponent of the density of states interpolates 
to 1 which is the Dirac limit. This limit is also obtained 
in the zero field limit in the isotropic case. Except for 
the zero energy finite value the energy dependence of the 
density of states in the case with disorder is similar to 
the lattice case. This suggests that the vortex disorder 
does not dramatically affect the density of states at low 
energies. An analysis of the IPR and the LDOS shows 
that the lowest lying states are delocalized, even though 
strongly peaked at the vortex cores. In the gapped s- 
wave case however the disorder introduces states in the 
gap thereby changing qualitatively the low energy den- 
sity of states, as in the high field limfljiS. We found a 
power law behavior with an exponent that scales with 
the magnetic length. 
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